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We investigate the exclusive electroproduction of ir + mesons from nucleons. To lead- 
ing twist, leading order a s accuracy the corresponding production amplitude can be 
decomposed into pseudovector and pseudoscalar parts. Both can be expressed in 
terms of quark double distribution functions of the nucleon. While the pseudovector 
contribution is connected to ordinary polarized quark distributions, the pseudoscalar 
part can be related to one pion t-channel exchange. We observe that the pseudovec- 
tor part of the production amplitude is important at xbj < 0.1. On the other hand, 
for 0.1 < XBj < 0.4 contributions from one pion exchange dominate. 
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1 Introduction 



Hard exclusive electroproduction of mesons from nucleons has become a field of growing 
interest. The recent progress in the QCD analysis of such processes is based on a factoriza- 
tion theorem |J (see also 0, || f|, |fl) valid in the case of longitudinally polarized photons: 
at large photon virtualities, Q 2 3> Aq CD , the underlying photon-parton sub-processes are 
dominated by short distances and, hence, can be calculated perturbatively. On the other 
hand, all information on the long distance dynamics of quarks and gluons can be collected 
in the distribution amplitude of the produced meson and in generalized (skewed |2|, || or 
double @) parton distribution functions of the nucleon. As a consequence, experimental 
investigations of meson production processes allow to probe details of the quark-gluon 
dynamics in nucleons beyond our current knowledge obtained mainly from inclusive high- 
energy processes. Ongoing and planned measurements at DESY (HERA M, HERMES 



i), CERN (COMPASS @), and Jefferson Lab [0 are therefore of great interest. 

In our previous work 0, |ll|] we have derived the leading order QCD production am- 
plitudes for neutral pseudoscalar mesons and vector mesons. The corresponding cross 
sections have been discussed within a specific model [E| for double parton distributions. 
In this paper we consider hard electroproduction of charged pions. Interest in this pro- 
cess arises in particular from the fact that it receives significant contributions from two 
production mechanisms: either (i) the photon interacts with a quark from the nucleon 
which, after hard gluon exchange, combines with a second quark of the target to the final 
pion, or (ii) the pion is produced from the meson cloud of the nucleon. The latter mech- 
anism is often referred to as the pion pole contribution. In meson-cloud models of the 
nucleon the second mechanism directly involves the electromagnetic form factor F n (Q 2 ) 



of the pion (see e.g. Jl3[ and references therein). Pion electroproduction has therefore 



often been considered as a generic process to determine F n (Q 2 ) (see e.g. |L4]]). However, 



as discussed in ref. |15[], there may be significant uncertainties due to type (i) contribu- 
tions. The latter have been modeled in ref. [15] in terms of hard gluon exchange diagrams 
with nonperturbative factors estimated by an overlap of light-cone wave functions with 
Chernyak-Zhitnitsky type distribution amplitudes [[17] . 

In this work we address questions similar to those raised in []15|| . The distinctive 
feature of our approach is a systematic use of perturbative QCD (PQCD) factorization 
0, |2|. Within a consistent PQCD analysis all necessary information about the nucleon 
structure is contained in quark double distributions. The latter are constrained through 
the behavior of ordinary quark distribution functions, as well as through sum rules and 
symmetry properties (for a discussion see ||, |16[). In particular we show that actually 
both production mechanisms, (i) and (ii), can be formulated in terms of quark double 
distribution functions of the nucleon, i.e. the existence of the pion pole contribution fits 
into the PQCD factorization framework. 

After modeling the involved double distribution functions we calculate the differential 
7r + production cross section in the region of small momentum transfers, —t < 0.5 GeV 2 . 
We find that both mechanisms, (i) and (ii), contribute significantly. While one pion 
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exchange dominates at intermediate values of Bjorken-x, perturbative quark exchange is 
relevant mainly at small x^y 

The paper is organized as follows: in Sec. 2 we present the amplitudes for charged 
pion electroproduction. Pseudovector and pseudoscalar parts of the relevant double dis- 
tribution functions are associated in Sec. 3 and 4 with the perturbative and pion pole 
contribution, respectively. Results are presented in Sec. 5. Finally we summarize. 



2 Amplitudes 

The calculation of meson electroproduction amplitudes is based on PQCD factorization 
[T[]. It applies to incident longitudinally polarized photons with large space-like momenta, 
Q 2 = —q 2 ^> Aq CD , and moderate momentum transfer to the nucleon target, \t\ < Aq CD . 
Using the techniques outlined in refs.|l|, 0, ^, |Tl| gives the following result for the ir + 
virtual photoproduction amplitude in leading twist, leading order a s accuracy: 
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Here the notation / d [x, y) . . . — f dx J£ dy . . . , with x = 1 — x is used. N(P, S) and 
N(P', S') are the Dirac spinors of the initial and scattered nucleon, respectively, with the 
corresponding four-momenta P, P' and spins S, S 1 . M stands for the nucleon mass. The 
average nucleon momentum is denoted by P = (P + P')/2, and the momentum transfer 
is r = P — P' with t = r 2 . The produced meson carries the four-momentum q' and 
q = (q + q')/2. Furthermore, we have introduced the variable to = 2q ■ P / (—q 2 )- Finally, 
n is a light-like vector with n ■ a = a + = a + a 3 for any vector a, and n = 7^71^. 

The long distance dynamics of the pion in the final state is contained in the decay 
constant f n = 133 MeV and the distribution amplitude $ w (r) (Ti 



(ir+(q')\Mx) l5 nMy)\0) = -ifAq' -n) C ^(r) e^^+^ dr , (2) 







where quark fields tp with proper flavor quantum numbers enter. 

The nucleon part of the amplitude (P is determined by double distribution functions 
AF du and AK du which are nondiagonal in flavor. They are defined as matrix elements of 



2 



a nonlocal quark operator sandwiched between proton and neutron states: 

(n(P', S>) | &(0) 75 n [0, z] M*) \p( p > 5 )>^=o = 

N(P', S') i 5 hN{P, S) J d [x, y] AF du + e «( p - z )-*(")AF rf j 

- N(P', S') 75 N(P, S) T -^ J d [x, y] ( e -**(^-*)-i»(*-*) AK du + j*{P-*)-*<r-*) AK du ) . 

(3) 

Here, as in the dependence of the double distribution functions on x, y and t has been 
suppressed. The up- and down-quark fields in (||) are separated by a light-like distance 
z ~ n. Gauge invariance is guaranteed by the path-ordered exponential 

[0,2] = Vexpl-igz,, [ A^zXjdX] 
Jo 

which reduces to 1 in axial gauge n ■ A = (g stands for the strong coupling constant and 
A 11 denotes the gluon field). 



3 Pseudovector contribution 



The production amplitude A w + contains a pseudovector and pseudoscalar contribution 
proportional to N(P', S') j 5 n N(P, S) and N(P', S') 75 N(P, S), respectively. We rewrite 
the pseudovector part according to isospin relations which connect flavor non-diagonal 



double distribution functions AF du with flavor diagonal ones. In particular we use |TT 

(n\ O du {z) \p) = (p\ 6 uu (z) \p) - (p\ 6 M (z) \p) , 



with 



6^'(z)=^ q {0)^n[0,z}^(z) 



z 2 =0 



(4) 



(5) 



In terms of flavor diagonal polarized double distribution functions the pseudovector 
part of the amplitude (^]) reads: 
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(6) 



Here e u and denote the electromagnetic charges of the up- and down-quarks. Note 
that Ak+ has a structure similar to the p + production amplitude derived in Oj: up to 
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prefactors the p + amplitude can be obtained from (|6|) by replacing polarized quark double 
distribution functions by unpolarized ones, with proper account of their opposite charge 
conjugation properties. 

In order to obtain a numerical estimate for the pseudovector contribution (|D to 7r + 
production models for the involved double distributions AF have to be constructed. We 
are guided here by the appropriate forward limit of double distributions, their symmetry 
properties, and sum rules which relate them to nucleon form factors. Following refs. [Z0| 
and we use: 



A: AF(x,y;t) = h(x, y) Aq(x) f{t), 

t y{l-x-y) 



B : AF(x,y;t) = h(x, y) Aq(x) exp 



A 2 x(l - x) 



(7) 



Here h(x,y) = 6y (1 — x — y)/{l — x) 3 , and Aq(x) denotes the corresponding ordinary 
polarized quark distribution. In numerical calculations presented in this paper we have 



used Gehrman- Stirling LO set-A parametrizations of polarized quark distributions [^1 
Furthermore, we have used the one loop expression for the running coupling constant with 
Aqcd — 200 MeV. Note that in addition to the above mentioned constraints, the models 
(|7|) are consistent with asymptotic solutions of evolution equations for double distributions 
and satisfy the symmetries of the latter Q. The t-dependence of double distributions in 
model A is governed by the form factor f(t) which, like in is taken to be equal to the 



nucleon pseudovector form factor [22 



The exponential dependence of AF on the momentum transfer t in model B can be 
obtained from a simple field theoretical investigation of double distribution functions 
outlined in [ffBl EIJ. For both models the scale A has been fixed at 1 GeV. With this 
choice model B reproduces the nucleon axial form factor up to momentum transfers 
t ~ —0.5 GeV 2 . Another important feature of model B is, that the small-x behavior of 
AF becomes less and less singular as \t\ increases. Such a behavior is expected since, 
like in the case of nucleon form factors, large momentum transfers filter out the minimal 
contribution to the Fock space wave function of the target. 



4 Pseudoscalar contribution 

As explained in ref.[^(J the double distributions AF(x, y\ t) can be related in the forward 
limit, t — 0, to ordinary polarized quark densities Aq(x). Therefore, it is possible to 
construct plausible models for AF, although little is known about their shape from first 
principles. On the other hand, for the double distributions AK no experimental infor- 
mation on corresponding forward densities is available which makes modeling in this case 
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more difficult. Nevertheless, some intuition about their magnitude has been obtained 



from model calculations p3[ . Still, contributions of the so-called .fT-terms to hard exclu- 
sive meson production have been neglected so far. In the present case of 7r + production 
the situation is different. As we will show, the pseudoscalar piece of the amplitude A n + 
can be associated with the virtual photoproduction of pions from the nucleon pion cloud. 
Due to the small pion mass this mechanism is expected to be important, especially at 
small momentum transfers t. 

We are interested in evaluating the contribution to the nucleon matrix element (|3|) 
arising from a situation when a quark-antiquark pair with low invariant mass propagates 
in the t-channel, such that it becomes close to a tt + bound state. To construct the corre- 
sponding double distribution function we introduce an effective Lagrangian describing a 
pseudovector pion- nucleon interaction [p2] : 

AriViv = ^ (0(x) 7 , 75 f0(x)) • (d^(x)) . (9) 

The nucleon and pion fields are denoted by <fi and (p, respectively, f are the common 
isospin matrices, and q^nn is the pion-nucleon coupling constant. Evaluating the matrix 
element (|3[) to first order in g n NN leads to: 

(n(P', S') | M0h 5 n [0, z] ^ u {z) \p(P, S)) z2=0 = V2g nNN N(P', S') 75 N(P, S) 

x J d A xe~ tr - x (O|T[^(O) 75 n[O,z]^^)0 + (x)] |0) 22=0 , (10) 

where ip + = ^75(^1+2 ¥2) is the pion field associated with the 7r + meson. The leading- twist 
contribution to the matrix element (|l(]) arises from the region of phase space where the 
quark fields are localized near the light-cone, z + ~ 0. Inserting a full set of intermediate 
pion states between the pion field and the non-local quark operator gives: 



(n(P', S')\ Vd(0)7 5 n [0, z\ ijj u {z) \p(P, S)) z2=0 

-i^/2g nN N 
ml — t 



N(P', S') l5 N(P, S) \ y *™ (0|^(0) 75 n [0, z] ^ u {z)\ir+{r)) z , =0 . (11) 



Using the pion distribution amplitude (§) and comparing with the definition of the pseu- 
doscalar double distribution functions in eq.(|J) yields: 

(AK d " + AK^)(x,y,tnO) = _ 2 ^ m 9,nn 5(x)%{yl (12) 
v ' — t 

As this contribution can be expressed entirely in terms of AK there is no danger of double 
counting if we add this term to the contribution arising from a model for AF, as discussed 
in the previous section. 

One can express the t-channel one pion exchange contribution also in the form: 

X + = -V2g nNN ^^P^^ Bl ■ W + r) F^+iQ 2 ), (13) 

m£ — t 
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where 

ei = ± (q' + r) (14) 
stands for the polarization vector of the longitudinally polarized photon, and 

is the leading QCD contribution to the n + electromagnetic form factor at large Q 2 |19| . 

To twist-2 accuracy, i.e. neglecting terms of order t/Q 2 and m 2 /Q 2 , the amplitude 
fll3|) is explicitly gauge invariant. This is guaranteed by the factor el ■ (q' + r) which 
arises from the hard photon-quark interaction. In the one pion exchange approximation 
the difference AK du — AK du decouples from the considered production process. Finally, 
note that in the case of ir° production a t-channel exchange contribution similar to fll2| ) 
appears in the production amplitude with a zero coefficient from the difference of two 
equal quark charges. 

To this end some comments regarding one pion t-channel exchange, as used above, 
are in order. In general one expects that at large center-of-mass energies any t-channel 
exchange of an hadronic state should be replaced by the corresponding Regge trajectory. 
Therefore, at small values of Bjorken-x additional contributions to the pseudoscalar part of 
the production amplitude A n + could occur. Once the pion pole contribution is replaced by 



reggeized pion (and rho) exchange p4[ the relation to the pion electromagnetic form factor 
seems to be lost (see, however, ref. [25] for an attempt to interpret pion electroproduction 
data in the framework of a Regge exchange model). We restrict ourselves in this work 
to the simplest model for the pseudoscalar part of the production amplitude in the hope 
that it contains the dominant contribution at small momentum transfers, i.e. close to the 
pion pole. 



5 Results 

In Figs.l and 2 we present results for the cross section of exclusive tt + production in the 
scattering of longitudinally polarized photons from nucleons at Q 2 = 10 GeV 2 . Up to 
logarithmic corrections the cross section is proportional to 1/Q 6 as one can infer from 
eq.®. 

In Fig.l we show the differential production cross section for the minimal kinematically 
allowed value of t = t min = — x|jM 2 /(l — £Bj)- We find that the pseudovector contribution 
(0) dominates at small values of the Bjorken variable XBj- Its XBj-dependence reflects the 
behavior of the relevant forward parton distributions ([?]). At x ~ 0.3 the cross section is 
governed by the pseudoscalar pion cloud contribution (p|). 

While models A and B give qualitatively similar cross sections at small momentum 
transfers, they start to differ for increasing \t\. This is illustrated in Fig. 2 where we 
present results for the differential production cross section taken at t = —0.4 GeV 2 . The 
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maximum value of Bjorken-x is determined by the condition \t\ > |t m i n |- For model A the 
pseudovector contribution decreases slowly with decreasing x B j and dominates the cross 
section at xbj ^0.1. At XBj ~ ^max the total pion production cross section is about half as 
large as obtained from the pion cloud contribution alone due to strong interference effects. 
On the other hand, in model B the a^-dependence of the pseudovector contribution is 
similar to the one from the pion cloud. Also, contrary to the previous case, at x B j ~ x max 
the production cross section is about 20% larger than obtained from the pion cloud alone. 

Finally we compare in Fig. 3 to the leading twist production cross section for ir° mesons 
which has been investigated in (J§. At t — t m[n the pseudovector contribution to ir + 
production is about a factor 5-10 larger as compared to the tc° case. At xbj ~ 0.3 large 
contributions from the pion cloud, which are absent in ir° production, add further to the 
cross section ratio. 

It should be emphasized that the dominance of the hard gluon exchange processes 
considered in this paper can only be ensured for sufficiently large values of Q 2 . Just like in 
the case of the pion form factor, competing contributions are provided by soft mechanisms 
which, in this case, correspond to the nonperturbative overlap of three hadronic wave 
functions. The interplay of soft and hard interaction mechanisms has been studied recently 



for the pion form factor in a QCD sum rule related model |pq| . Guided by the results of 
this investigation we expect the dominance of hard gluon exchange for values of Q 2 above 
10 GeV 2 . 



6 Summary 

We have discussed the hard exclusive electroproduction of ir + mesons from nucleons. For 
longitudinally polarized photons one can express the corresponding photoproduction am- 
plitude to leading twist, leading order a s accuracy in terms of quark double distribution 
functions of the nucleon. We find that the amplitude receives contributions from pseu- 
dovector and pseudoscalar pieces. The pseudovector contribution is determined by quark 
double distributions with well known properties in the forward limit, as given by ordinary 
polarized quark distribution functions. However, for the pseudoscalar part such informa- 
tion is not available. In this work it has been modeled in terms of an interaction of the 
virtual photon with the nucleon pion cloud. As a consequence, the related double distri- 
bution is determined by the pion distribution amplitude. The corresponding contribution 
to the pion production amplitude turns out to be proportional to the pion form factor. 

We have found that both, pseudovector and pseudoscalar terms should be included 
if one considers the whole < XBj < 1 range. In particular, we have observed that the 
pseudovector part dominates at small XBj < 0.1. On the other hand, the pion cloud 
contribution controls the n + production at small t ~ t min and 0.1 < XBj < 0.4. At 
larger values of \t\ ~ 0.4 GeV 2 the relative weight of pseudoscalar and pseudovector 
contributions has been found to be sensitive to details of the relevant double distribution 
functions. Further investigations of double distribution functions, as well as experimental 
data on exclusive meson production are certainly needed. 



7 



Acknowledgments: 

G.P. would like to thank the SLAC theory group for their hospitality. Furthermore, we 
would like to acknowledge discussions with V. Burkert, S.J. Brodsky, L. Frankfurt, M.V. 
Polyakov and P. Stoler. The work of A.R. was supported by the US Department of Energy 
under contract DE-AC05-84ER40150. 

After this paper has been completed, we learned about an independent calculation of 
7r + electroproduction which uses skewed parton distributions computed in a chiral quark 
soliton model |27J]. The corresponding results turn out to be quite similar to ours. We 
thank the authors of [^7[ for discussions and an exchange of results. 



References 

J.C. Collins, L.L. Frankfurt, M. Strikman, Phys. Rev. D 56, 2982 (1997). 

A. Radyushkin, Phys. Lett. B 380, 417 (1996); Phys. Lett. B 385, 333 (1996); Phys. 
Rev. D 56, 5524 (1997). 

P. Hoodbhoy, Phys. Rev. D 56, 388 (1997). 

L. Mankiewicz, G. Piller, T. Weigl, Eur. J. Phys. C 5, 119 (1998). 
M. Vanderhaeghen, P. A.M. Guichon, M. Guidal, Phys. Rev. Lett. 80, 5064 (1998). 
X. Ji, J. Phys. G 24, 1181 (1998). 

JA. Crittenden, Exclusive production of neutral vector mesons at the electron- 
positron collider HERA, Springer, Berlin 1997 |[hep-ex/970 100S ] . 



G. van der Steenhoven, Future physics with light and heavy nuclei at HERMES, in 
Proc. Workshop on Future Physics at HERA, ed. G. Ingelman et ai, DESY, Hamburg 
1996. 

G. Baum et al, COMPASS proposal, Report No. CERN-SPSLC-96-14. 
See e.g. C. Carlson, Nucl. Phys. A 622, 66c (1997). 

L. Mankiewicz, G. Piller, T. Weigl. Hard leptoproduction of charged vector mesons, 
Report No. TUM-T39-98-12, hep-ph/9712508. 

A.V. Radyushkin, Phys. Rev. D 59, 014030 (1998). 

J. Speth and V.R. Zoller, Phys. Lett. B 351, 533 (1995). 

C.J. Bebek et al, Phys. Rev. D 17, 1693 (1978). 

C.E. Carlson and J. Milana, Phys. Rev. Lett. 65, 1717 (1990). 



8 



[16] A.V. Radyushkin, Symmetries and structure of skewed double distributions, Report 
No. JLAB-THY-98-41, hep-ph/9810466. 

[17] V.L. Chernyak and A.R. Zhitnitsky, Phys. Rep. 112, 173, 1984. 

[18] A.V. Efremov and A.V. Radyushkin, Phys. Lett. B 94, 245 (1980). 

[19] S.J. Brodsky and G.P Lepage, Phys. Lett. B 87, 379 (1979). 

[20] A.V. Radyushkin, Phys. Rev. D 58, 114008 (1998). 

[21] T. Gehrmann and W.J. Stirling, Phys. Rev. D 53, 6100 (1996). 

[22] T. Ericson and W. Weise, Pions and Nuclei, Oxford University Press, 1988. 

[23] X. Ji, W. Melnitchouk, X. Song, Phys. Rev. D 56, 5511 (1997); 

V. Petrov, P. Pobylitsa, M. Polyakov, I. Bornig, K. Goeke, C. Weiss, Phys. Rev. D 
57, 4325 (1988). 

[24] M. Guidal, J.-M. Laget, M. Vanderhaegen, Nucl. Phys. A 627, 645 (1997). 

[25] M. Vanderhaegen, M. Guidal, J.-M, Laget. Phys. Rev. C 57, 1454 (1998). 

[26] A. Szczepaniak, C.-R. Ji, A.V. Radyushkin, Phys. Rev. D 57, 2813 (1998). 

[27] K. Goeke, M. Penttinen, M.V. Polyakov, RUB-TPII-19/98; 

L.L. Frankfurt, PV. Pobylitsa, M.V. Polyakov, M. Strikman, RUB-TPII-20/98. 



9 




Figure 1: Differential production cross section for exclusive tt + production through the 
scattering of longitudinally polarized photons from protons at t — t min and Q 2 = 10 GeV 2 . 
The dotted and dot-dashed curves show the pseudovector and pseudoscalar contributions, 
respectively. The upper figure corresponds to model A, the lower one to model B. 
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Figure 2: Differential production cross section for exclusive tt + production through the 
scattering of longitudinally polarized photons from protons at t — —0.4 GeV 2 and Q 2 = 
10 GeV 2 . The dotted and dot-dashed curves show the pseudovector and pseudoscalar 
contributions, respectively. The upper figure corresponds to model A, the lower one to 
model B. 
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Figure 3: Ratio of the differential production cross sections for exclusive ir + and 7r° 
production through the scattering of longitudinally polarized photons from protons at t — 
t min and Q 2 = 10 GeV 2 for model A. The dotted curve corresponds to the pseudovector 
contribution alone. 
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